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Abstract 

Using dynamical mean-field theory, we have evaluated the magnetic instabilities and T = phase 
diagram of the double-exchange model on a Bethe lattice in infinite dimensions. In addition to 
ferromagnetic (FM) and antiferromagnetic (AF) phases, we also study a class of disordered phases 
with magnetic short-range order (SRO). In the weak-coupling limit, a SRO phase has a higher 
transition temperature than the AF phase for all fillings p below 1 and can even have a higher 
transition temperature than the FM phase. At T = and for small Hund's coupling Jh, a SRO 
state has lower energy than either the FM or AF phases for 0.26 < p < 1. Phase separation is 
absent in the Jh — > limit but appears for any non-zero value of Jh- 
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Due to the keen interest in manganites [| and dilute magnetic semiconductors the 
double-exchange (DE) model remains one of the central models in condensed-matter science. 
Yet the behavior of the DE model for small Hund's coupling continues to puzzle researchers. 
Both the nature of the uniform phases and the presence of "phase separation (PS) in the 
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weak-coupling regime is the subject of ongoing debate |3, |4|, Hj |6( . This paper uses dynamical 
mean-field theory (DMFT) to obtain the magnetic instabilities and T = phase diagram of 
the DE model in infinite dimensions. We reach the surprising conclusion that a short-range 
ordered (SRO) phase is stable for small Hund's coupling J H and even has a higher transition 
temperature than the long-range ordered ferromagnetic (FM) and antiferromagnetic (AF) 
phases for some fillings. 

Developed in the late 1980's by Muller-Hartmann []| and Metzner and Vollhardt 0, 
DMFT exploits the momentum independence of the self-energy in infinite dimensions. Even 
in three dimensions, DMFT is believed to capture the physics of correlated systems including 



the narrowing of electronic bands and the Mott-Hu 



Dard 



tranati<m,|gL Although DMFT 
until now there 



has been widely applied to the DE model fl J fi, Q lj Q HQ, 
has been no complete treatment of the magnetic instabilities and T = phase diagram of 
the DE model. Usually, DMFT calculations are performed on either a Bethe or hypercubic 
lattice. Because of the simplified formalism on the Bethe lattice and the risk of pathological 
results jljj associated with the unbounded, hypercubic density-of-states (DOS), we choose 
to work on a Bethe lattice with a semi-circular DOS bounded by ±W/2. 
The DE Hamiltonian is 



H = -tJ2 (4a C ja + c} a C ia ) - 2 J H Yl S i ' S * i 1 ) 

where c\ a and Ci a are the creation and destruction operators for an electron with spin a 
at site i, = (l/2)c\ a a al 3C i/3 is the electronic spin, and Sj = Snij is the spin of the local 
moment (treated classically). Repeated spin indices are summed. Then within DMFT, the 
effective action on site is given by 

Scff(m) = -T^co a (zV n ){<3 (zV n );^j + Ja af3 ■ m}c 0( g(wn), (2) 

n 

where J = JrS, v n = {2n + l)irT, Co a (iu n ) and CQ a {iv n ) are now anticommuting Grassman 
variables, and G (iu n ) is the bare Green's function containing dynamical information about 
the hopping of electrons from other sites onto site 0. 



In the infinite-dimensional limit, the high-temperature non-magnetic (NM) phases of the 
Heisenberg and DE models have a vanishing correlation length £. The SRO phase is a 
bulk solution of the DE model that, like the FM and AF states, relies only on the local 
topology of the Bethe lattice. It has some of the same characteristics as conventional spin 
glasses: a finite local magnetization and spin-spin correlations that decay exponentially over 
distance The SRO phase is characterized by a correlation parameter q that gives the 
average q = (sin 2 (^/2)), where Q{ is the angle between a central spin and each neighboring 
spin. Overall, the neighboring spins describe a cone with angle 2arcsin( v /g) around the 
central spin. The FM and AF phases have, respectively, q = and 1. If M n is the average 
magnetization of the lattice sites at a distance na (a is the lattice constant) from the central 
site, then M n = (1— 2g)M n _i so that the magnetization about every site decays exponentially 
like \M n \ = \Mq\ exp(— na/£). The correlation length £ = — a/log|2g— 1| diverges only in 
the FM and AF limits and vanishes in the NM state obtained by setting q = 1/2. 

In their earlier work on the Bethe lattice, Chattopadhyay et al. |5( assumed that the 
angles Qi were the same for every neighbor and inaccurately characterized this state by its 
"incommensurate correlations" rather than by its short-range order. Since only the FM 
and AF states possess well-defined wavevectors on the Bethe lattice, the SRO state cannot 
be interpreted using the formalism originally developed by Muller-Hartmann |7] for the 
hypercubic lattice 9(. Ordered spiral solutions are not allowed on the Bethe lattice due to 
its topology and lack of translational invariance. 

Because the SRO phase has no well-defined wavevector, its transition temperature cannot 



be obtained from the magnetic susceptibility 
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Rather, Tsro(<?) must be solved from 



coupled Green's function equations. On a Bethe lattice, the bare and full local Green's 
functions are related by 

W 2 

G^(z^)- 1 = z n L- —{qGW(iu n ) + (1 - q)GM(iv n )}, (3) 

where G {v) and are spin-reversed Green's functions and z n = iu n + /i. The full Green's 
function Gj v '(ii/ n ) at site is obtained using the effective action of Eq.(J2J). With the param- 
eterization 

G$\iv n y l = (z n + iQl+Q^Gz, (4) 

we evaluate the critical transition temperatures by linearizing in the magnetization M = 
(m z ) on site (taken to be spin up). To first order, R n is independent of M and is obtained 




FIG. 1: The correlation parameter q (solid) and the associated maximum SRO transition temper- 
ature 7g^Q X ^ (dashed) versus filling p for a SRO state in the weak-coupling limit. Also plotted in 
the thin short-dashed line are the weak-coupling limits of Tq for q = and Tn for q = 1. 

by solving the cubic equation 




while Q^P = —Q$ is proportional to M. After integrating exp(— (m)) over the Grassman 
variables, we find that the probability for the local moment to point along m is P(m) oc 
exp((3J eff Mm z ) with J eff = -2T J J2 n {Q^ / M) / {{z n + R n ) 2 - J 2 ). Hence, T SRO {p,q) is 
solved from the implicit relation Tsro = </ e //(7sRo)/3 or 

2 J 2 R 
1 = ~~ (2<? ~ X) ? (z n + R n ) 2 (z n + 2(1 - q)R n ) - J 2 (z n + 2(2 - q)R n /3) ' (6) 

which correctly reduces to the FM or AF results when q = or 1. Notice that TsrqQo, q) = 
in the NM state with q = 1/2. Of course, Tsro(p, q) is the same for any arrangement of 
the spins with the same average value of sin 2 (#j/2) . In terms of the chemical potential 
/i, the electron filling p is given by p = 1 — (32T/W 2 ) J2 n R- e ^n- So p = 1 when /i = 0, 
corresponding to half filling or one electron per site. 

In the weak-coupling limit, T <^ J W, TsRo(p,q) is evaluated by converting the 
Matsubara sum in Eq.© into an integral over v and evaluating R(z = %v + fi) to zeroth 
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order in Jh- Maximizing Tsro(p, q) with respect to q, we obtain TgRQ (p) plotted in Fig.l. 
Due to the symmetry about p = 1, we restrict consideration to values of p between and 
1. Remarkably, a SRO phase with AF correlations (1/2 < q < 1) can always be found with 
a higher transition temperature than the AF phase. As q — *■ 1, the AF and SRO transition 
temperatures meet at p = 1 or half filling. Even in the range of fillings between about 
0.26 and 0.38, where the Curie temperature is nonzero, a SRO phase with FM correlations 
(0 < q < 1/2) has the higher critical temperature! An instability is also found in the range 
of fillings between 0.38 and 0.73, where neither the FM nor AF phases are stable for small 
J/W. The point at which q — 1/2 and T^^\p) = lies slightly below p = 0.5. The 
divergence of the weak-coupling results as p — > 1 signals the breakdown of the condition 
T ^ J under which they were derived and is associated with the appearance of a gap in 
the AF DOS, as discussed below. The correlation parameter q changes discontinuously at 
TgRO (p) from 1/2 in the NM phase above to a value less than or greater than 1/2 in the 
SRO phase below. 

We have also obtained the magnetic instabilities of the DE model for arbitrary J/W. 
For any J/W and p, Fig. 2 (a) indicates the phase with the highest transition temperature. 
Around p = 1, an AF instability occurs in a narrow range of fillings that vanishes as 
J/W —>■ or oo. The range of AF instabilities has a maximal extent (from p = 0.78 to 1) 
when J/W « 0.33. For J/W < 0.33, T^o\p) = when q = l/2 along the NM curve. For 
0.33 < J/W < 0.5, Tg^Q (p) > for all fillings and q jumps from in the FM to a value in 
the range 1/2 < q < 1 in the SRO phase. For J/W > 0.5, either the FM or AF phase has a 
higher transition temperature than the SRO phase and q jumps from to 1 at the FM/AF 
boundary. 

The ground-state energies of the FM, AF, and SRO phases are obtained after once again 
parameterizing Gq , by Eq.pjl. Then R(z) solves the quartic equation 

W 2 

R(z + R){(z + 2(1 - q)Rf - J 2 } + — {z + 2(1 - q)Rf = 0, (7) 

which differs from Eq.(J5|) (except when q — 1/2) because we are now working in the broken- 
symmetry phase. The interacting DOS per spin is defined by 

which reduces to the bare DOS iVo(e) = (8/ttW 2 ) ^W 2 /A - e 2 when J -> 0. The 
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FIG. 2: The (a) magnetic instabilities and (b) T = phase diagram of the DE model. 

AF DOS with q — 1 can be solved analytically for all J: N(u, q = 1) = 
(8\uj\/nW 2 )Re ^J{W 2 /4: + J 2 - uj 2 )/(uj 2 - J 2 ), which vanishes for \u\ < J and \u\ > 
\jw 2 /4 + J 2 . Hence, N(u, q = 1) contains a square-root singularity on either side of a 
gap with magnitude 2J. As J/VF — > oo, the width of each side-band narrows like W 2 /8J. 
It is also possible to explicitly evaluate the DOS at u = for any q: N(ou = 0,q) = 
(8/nW 2 (l - q))Re^/j 2 - J 2 , which vanishes when J > J c = W(l - q)/2. In a NM, the 
critical value required to split the band is J c = W/4, as found earlier In an AF, J c = 
since a gap forms for any nonzero coupling constant. 

Taking J/W = 0.2, we plot the interacting DOS versus u for q = 0,1/4,1/2,3/4, and 
1 in Fig.3. The FM DOS has kinks at uj = ±(W/2 — J). Between the kinks, both up- 
and down-spin states appear; on either side, the bands are fully spin-polarized. The kinks 
disappear in the SRO phase with 1 > q > due to the absence of long-range magnetic 
order. As q increases from to 1, the width of the DOS shrinks and a gap appears when 
J > J c (q)- The breakdown of the weak-coupling results in the limit p — > 1 is caused by the 
appearance of a gap in the AF DOS for any nonzero J. 
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FIG. 3: The interacting DOS N(u,q) (normalized by 1/W) versus tu/W for J/W = 0.2 in the 
FM or q = (solid), SRO phase with q = 1/4 (long dash), NM or q = 1/2 (dash), SRO phase with 
q = 3/4 (small dash) and AF or q = 1 (dot dash) phases. 

Quite generally, the total energy E(p, q) of any state may be written as an integral over 
the interacting DOS: E(p,q)/N = 2 J du uo f(uo)N(uo,q), where f(u) = l/(exp(f3(u -//)) + 1) 
is the Fermi function and \x = n{p). At T = 0, this relation becomes 

1 , . If, { 16Rz) 

^,g) = --y^Re|l + — j, (9) 

which converges because — > — W 2 /(16z) as |d — * oo. 

In their numerical work, Chattopadhyay et a/. [5] obtained a very complex phase diagram 
at T = 0, with SRO phases above J/W ~ 1/16 and PS between AF and FM phases for 
0.1 < p < 1 below this value. In light of the numerical difficulty of constructing the phase 
diagram for small J/W, we have studied this limit analytically. As demonstrated below, PS 
disappears in the limit J/W — > 0. 

Carefully accounting for the dependence of the chemical potential fi(p) on J/W for a 
: ixed filling p, we find that the difference between the energies of the SRO and NM phases 
16^ is AE(p, q)/N = — 3T SRO (p, <?)/2, where both sides are evaluated analytically to second 
order in J/W. So for small J/W, the ground state energy is minimized by the same corre- 
lation parameter q that maximizes the transition temperature! This result is not surprising: 
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for small J/W, the transition temperature is also small so the correlation parameter that 
appears at Tsro(p, q) continues to minimize the energy at T = 0. 

It is now easy to understand why a SRO phase has a lower energy than the AF phase 
for small J/W when p < 1. Due to the narrowing of the AF DOS, the AF energy may be 
higher than that of a SRO phase with J < J c (q) and no energy gap. When the chemical 
potential of the AF lies outside the energy gap, then the square-root singularity of the AF 
DOS will raise the energy compared to the broader, gapless DOS of a SRO state. 

For any filling p, the energy E^ mm " l (p) is obtained by minimizing E(p,q) with respect to 
q. When J > 0, PS occurs due to the formation of an energy gap in the AF phase, so 
that (l/N)dE^/dp\ p=1 - = - J. A necessary condition for PS between fillings p\ < 1 and 
P2 = 1 is that the second derivative of the energy E^ min \p) must change sign at some filling 
p* between p\ and pi- To second order in J/W, the condition d 2 E^ mm \p) / dp 2 = may 
be written l/\Jl — (2fi/W) 2 oc (J//i) 2 . Since the right-hand side vanishes as Jh —* 0, the 
condition for PS cannot be satisfied except as /i — > or p* — > 1 . Near p = 1, fx/W oc 1 — p 
so the PS width p2 — p\ grows linearly with J/W. 

The magnetic instabilities and T = phase diagram of the DE model are plotted side-by- 
side in Fig. 2. For J/W < 0.33, a SRO phase is stable over a range of fillings that diminishes 
with increasing J/W. The NM solid curve in the SRO region with q = 1/2 separates a 
SRO phase with AF correlations (1/2 < q < 1) on the left from one with FM correlations 
(0 < q < 1/2) on the right. For J/W < 0.33, PS occurs between a SRO state and an AF 
with p = 1. In agreement with the discussion above, the PS region grows linearly with J/W. 
For J/W > 0.33, the SRO phase is bypassed and PS is found between a FM and an AF 
with p = 1. The horizontal dashed line in Fig. 2(b) separates these two PS regions. The PS 
width diminishes as J/W increases past 0.33. In this high-coupling regime, our results are 
in qualitative agreement with other authors . Aside from the behavior at very small 

J/W, Fig. 2(b) also agrees well with the phase diagram in Ref.J^j]. 

Together, Figs. 2 (a) and (b) provide a bird's eye view of the magnetic instabilities and 
phase evolution in the DE model. As the temperature is lowered, the phase space of the pure 
FM shrinks and the AF phase with p < 1 disappears altogether. The coupling J/W ~ 0.33 
below which a SRO phase is stable at T = closely agrees with the value where the NM 
curve intersects the FM instability boundary and where the AF instability covers the widest 
range of fillings in Fig.2(a). Bear in mind that the pure FM, AF, and SRO phase instabilities 
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of Fig.2(a) may be bypassed when {J/W,p} lies in the PS region of Fig.2(b). 

The presence of a NM or SRO phase at T = would seem to contradict Nernst's theo- 
rem, which requires that the entropy is quenched at zero temperature. However, entropic 
ground states occur quite frequently in infinite dimensions, appearing in earlier work using 
a hypercubic lattice for the frustrated Hubbard l| and DE models. 

As pointed out over 40 years ago the competing FM and AF interactions in the 
DE model produce a RKKY-like interaction for small J/W jf]], which frustrates the ordered 
phases and stabilizes a SRO state. For small J/W, the component of the electronic spin 
along the local- moment direction is of order J/W so the electronic spins are not frozen at 
T = 0. When the dimension is lowered, the SRO state evolves into the state with incommen- 
surate correlations (IC) obtained by Monte Carlo simulations |3|. Indeed, simply replacing 
"SRO" by "IC," Fig. 2(b) bears a striking resemblance to the one- and two-dimensional phase 
diagrams of Ref . |3[ . 

Whether the SRO phase is a new kind of spin glass or spin liquid can only be resolved 
by future studies. Two of the most important unresolved questions about spin glasses [l^ 
are whether there exists a true thermodynamic transition and whether a model without 
quenched disorder can support a spin glass The answers for the SRO phase of the 

DE model on a Bethe lattice are clear. Although not marked by a divergent susceptibility, 
the SRO transition is characterized by the development of short-range magnetic order and 
a reduction in the entropy compared to the NM state. In future work, we will pursue the 
analogy with the spin-glass solution for the random Ising model on a Bethe lattice, where 
the Edwards- Anderson order parameter can be explicitly constructed j^ . 

To summarize, we have studied the behavior of a SRO solution to the DE model in infinite 
dimensions. Remarkably, the SRO transition temperature may be higher than that of the 
ordered FM and AF phases and the SRO phase remains stable for small couplings down to 
T = 0. The possibility of analytically constructing a disordered solution for the DE model 
in infinite dimensions should be a boon to future investigations of the DE model in lower 
dimensions. 
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